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Although  direct  measurements  of  these  dynamic  properties 
ABSTRACT  seem  feasible,* ’*  they  are  not  available  at  this  time. 

Thus  theoretical  calculation  of  these  values,  until  data 
A theoretical  study  of  mass  movement  on  the  Mis-  become  available,  must  be  invoked.  Also,  the  importance 
sissippi  River  delta  front  has  been  made  using  recently  ^ theoretical  model  as  a tool  to  aid  experimentation 
acquired  field  data  and  a simple  rheological  model,  should  not  be  neglected.  Models  define  a set  of  specific 
Recent  measurements  of  sediment  properties  include  cohe-  parameters  that  would  need  to  be  measured  and  the  accu- 
sion,  bulk  density,  pore  pressure,  and  internal  friction  racies  that  are  needed.  This  information  could  greatly 
angle.  Also,  the  geometry  of  a typical  type  of  instabil-  improve  the  efficiency  of  instrument  design,  data  col- 
ity  feature,  an  elongate  slide,  is  examined  in  detail.  A lection,  and  ^ata  analyses, 
rheological  model  describing  a Coulomb-viscous  plastic 

in  effective  stress  terms  is  proposed  to  describe  cer-  CHARACTERISTICS  OF  SUBMARINE  MASS  MOVEMENTS 
tain  features  of  the  mass  movement  process. 

A variety  of  slope  and  sediment  instability  mor- 
The  model  is  used  to  derive  equations  defining  the  phologies  have  been  identified  and  have  been  interpreted 
initial  failure  of  the  slopes,  the  mass  movement  thick-  collapse  depressions,  bottleneck  slides,  elongate 
ness  and  accelerations,  the  velocity  of  fl^  in  the  gUjes,  shallow  rotational  slumps,  and  mudflow  gullies 
gullies,  and  the  shape  of  the  mud  nose  scarp  in  the  toe  overlapping  depositional  lobes."*  A schematic  dia- 

area.^  The  model  indicates  the  importance  of  pore  pres-  gj-a,,  of  jhe  features,  taken  from  Roberts  et  al.,*  is 
sure  in  controlling  sediment  dynamics.  shown  in  Fig.  l!  Although  each  of  the  features  possesses 

somewhat  distinct  surface  morphologies,  certain  types 
Results  of  the  model  compare  favorably  with  the  app^g^  to  grade  into  others  or  to  be  linked  spatially, 
observed  shapes  of  mud  noses.  Initial  failures  on  low  fhis  suggests  that  there  may  be  a great  deal  in  common 
slopes  (~0.5®)  are  attributed  to  pore  pressures  ap-  among  the  physical  factors  that  cause  sediment  failure 
proaching  geostatic  values.  Flow  velocities  are  calcu-  and  movement, 
lated  to  be  several  feet  per  second,  based  upon  esti- 
mated sediment  viscosities  and  pore  pressures  during  Examination  of  the  various  types  of  bottom  fea- 

tures  indicates  that  the  most  common  and  extensively 
occurring  type  is  the  elongate  slide  (Fig.  2).  The  slide 
INTRODUCTION  three  basic  components;  (1)  a somewhat  circular 

area  of  multiple  scarps  and  blocks,  leading  into  (2),  a 
Recent  detailed  mapping  of  the  Mississippi  delta  long,  linear  gully,  and  terminating  in  (3),  a composite 
front  has  revealed  the  widespread  occurrence  of  bottom  I®”  lobe  area.  The  overall  length  of  the  slide 

features  that  appear  to  result  from  sediment  mass  move-  ™®y  mach  several  thousand  feet.  The  slump  areas  show 
ments  .***’* ’^  The  stability  of  these  sediments  has  been  curved  scarps  that  vary  from  10  to  25  ft  in  height  and 
generally  discussed  using  the  concept  of  effective  have  slopes  of  1®  to  A®.  The  basal  shear  planes  of  the 
stress;*’*’^  however,  there  have  been  no  attempts  to  scarps  are  only  slightly  concave  upward  and  tend  to 
analyze  the  dynamic  aspects  of  the  problem  in  this  merge  at  depth  into  a single  shear  surface  that  is 
light.  This  report  presents  the  results  of  a theoreti-  inclined  parallel  to  the  sediment  surface.  The  average 
cal  study  that  used  recently  acquired  field  data  and  a depth  of  the  basal  shear  plane  is  less  than  about  100  ft, 
simple  rheological  model  to  examine  the  dynamics  of  mass  and  its  slope  is  about  0.5®.  The  slide,  which  begins  as 
movements  on  the  delta  front.  a shallow  rotational  slump,  quickly  becomes  transla- 

tional, showing  substantial  downslope  displacements. 


Although  the  geometry  of  the  bottom  features  has 

been  described,  present  understanding  does  not  indicate  

the  role  of  pore  pressures  in  determining  how  movements  References  and  illustrations  at  end  of  paper, 
take  place,  what  the  speeds  or  depth  of  flow  may  be,  or 
how  sediment  properties  may  change  during  movement . 


determined  for  delta  sediments,  these  measurements  have 
not  been  coordinated  or  systematically  taken.  With  this 
lack  of  a proper  data  set,  the  application  of  the 
theoretical  results  to  submarine  landslides  must  be 
considered  exploratory. 

Initial  failure.  From  the  geometry  of  the  fea- 
tures that  have  been  identified,  an  initial  analysis  of 
the  mechanisms  of  these  subaqueous  slope  instabilities 
can  be  made  in  terms  of  equation  (S),  which  at  failure  is 

C - yZ  (1  + p/yZ)  ()> 

1 — s 9 « ^ 


It  is  possible  Co  calculate  from  this  equation  the  pore- 
water  pressure  (p)  needed  to  initiate  failure.  The 
general  slope  angle  for  the  delta-front  features  ranges 
from  0.1®  to  1.1®,  and  the  unit  weight  of  the  delta-front 
clays  is  approximately  100  Ib/ft*.^  Calculations  of 
pore-water  pressure  needed  for  failure  have  been  made 
using  various  values  of  C and  R and  for  a thickness  of 
sediment  of  SO  ft.  This  depth  represents  a reasonable 
estimate  of  the  depth  of  the  basal  shear  plane. 

Table  1 shows  calculated  values  of  pore-water 
pressure  needed  for  failure.  The  pressures  are 
expressed  as  the  constant  G or  the  pore-water  pressure 
ratios  to  geostatic  pressure  p/fZ.  For  all  cases,  the 
pore-water  pressure  needs  to  be  very  large  for  failure 
to  occur.  It  must  be  in  excess  of  hydrostatic  pressure, 
and  this  constitutes  strong  artesian  pressure.  In  addi- 
tion, it  is  clear  that  the  calculated  values  approach 
(or  very  slightly  exceed)  geostatic  pressures,  repre- 
senting a condition  of  almost  zero  effective  stress. 

Preliminary  results  of  pore-water  pressure  mea- 
surements within  the  subaqueous  sediments  on  the  Missis- 
sippi delta-front  slope  have  been  reported.*  The 
piezometers  were  installed  in  block  28,  South  Pass  area, 
at  depths  of  about  26  ft  and  49  ft  below  the  mudline. 
The  pore-water  pressure  data  revealed  large  excess  pres- 
sures after  7 hours  of  stabilization  of  the  system,  and 
values  were  considerably  in  excess  of  hydrostatic  pres- 
sure. Significantly,  the  ratios  of  p/yZ  approach  the 
geostatic  condition,  with  values  of  0.895-0.946.  Even 
larger  pressures  have  also  been  recorded  elsewhere  in 
the  offshore  delta  region,  with  p/yZ  ratios  of  0.986  at  a 
depth  of  50  ft  within  the  sediment  (Bennet,  personal 
communication) . 

Thus,  direct  evidence  indicates  that  large  pres- 
sures occur  within  the  delta-front  slope  sediments. 
Comparison  of  the  measured  pore-water  pressure  ratios 
(p/yZ)  with  the  calculated  ratios  needed  for  failure 
shows  very  close  agreement.  Indeed,  using  " 0,  • 

20®  (Henkel)  and  the  highest  ratios  actually  recorded, 
failure  can  be  achieved.  With  larger  values,  only 

relatively  small  increases  in  pore-water  pressure  are 
needed  for  failure.  For  example,  for  • 200  Ib/sq  ft 

and  B ■ 0.5®  the  pore  pressure  has  to  be  increased  by 
about  500  Ib/aq  ft  for  failure  to  occur.  Fluctuations  of 
150  to  300  Ib/sq  ft  have  been  recorded  resulting  from  the 
passage  of  waves  (Bennet,  personal  communication). 

The  spatial  distribution  of  the  landslides*  seems 
to  reflect  the  fact  that  primary  failures  are  generally 
localized  in  shallow  water  or  downslope  of  the  deposi- 
tional  lobes  of  landslides  in  shallower  water.  Further- 
more, the  pore  pressures  need  to  be  at  or  exceeding 
geostatic  pressures.  Pore  pressures  after  failure  may 


then  be  expected  to  increase  further,  exceeding  geo 
static  values. 

Flowage  of  sediment  in  gullies.  The  behavior  of 
sediments  in  the  chute  areas  can  also  be  addressed  using 
theoretical  results.  Given  that  the  excess  pore  pres 
sures  which  built  up  and  exceeded  the  failure  criteria 
are  maintained,  then  the  average  velocity  of  the  flows 
can  be  calculated.  Initially  sediment  speed  would  in- 
crease until  quasi-equilibrium  condition  was  reached. 

The  velocity  change  during  acceleration  can  be 
determined  from  equation  (11).  Taking  as  an  example  8 ■ 
0.5®,  C =•  lO.Olb/sq  ft,  h.  - 50  ft,  41  - 20®,  and  G - 

8 D 8 

1.0  , the  acceleration  has  a maximum  value  of  .216 

ft/sec*.  After  only  90  sec  of  acceleration  the  velocity 
reaches  a maximum  value  of  6.5  ft/sec.  These  velocities 
are  highly  dependent  upon  the  pore  pressures.  When  G “ 
,99,  the  maximum  velocity  is  2.4  ft/sec. 

The  parameter  e will  play  an  important  role  in 
determining  the  flow  characteristics  of  the  sediments. 
It  is  a function  of  the  bottom  slope,  friction  angle,  and 
pore  pressure  ratio  to  geostatic  pressure.  Thus 


0 - (1  - G) 


where  G “ ^ 

The  lower  limit  of  the  range  of  values  of  e is  determine 
by  the  failure  criterion.  From  Table  1,  the  ratio  G at 
failure  had  values  from  0.98  to  about  1.14,  for  B * 0.50 
and  4)  = 20®.  Therefore,  e has  minimum  values  ranging 
from  about  0.001  to  0.06.  Thus,  during  flow,  values 
greater  than  these  are  expected.  The  upper  limit  o'  e 
obtainable  is  not  easy  to  define;  however,  for  a pore 
pressure  ratio  to  geostatic  pressure  of  1.25  e has  the 
value  0.100.  Therefore,  the  range  of  variation  of  C may 
be  taken  as  between  about  0.0  and  O.I.  Once  the  slope 
and  the  friction  angle  41  are  fixed,  C is  directly 
proportional  to  the  pore  pressure  ratio  G. 

The  properties  of  the  sediment  flow  within  the 
chute  areas  can  be  estimated  from  equations  (8),  (9), 
and  (10).  Assume  that  the  chute  sediments  can  be  charac- 
terized by  8 " 0.5®,  *10,0  Ib/sq  ft,  4' ^ * 20®,  h^^  • 

50  ft,  and  y =•  100  Ib/cu  ft.  The  thickness  of  the  plug 
or  rafted  material,  from  equation  (8),  is  33.3  ft  for  t 
« .003  and  40  ft  for  e • .0025.  Thus  most  of  the 
sediment  column  within  the  chute  would  be  unsheared, 
with  the  shear  zone  concentrated  in  the  bottom  10  ft  to 
17  ft.  If  the  point  is  reached  where  e decreases  to  its 
minimum  value,  in  this  example  .002,  then  ■ 50  ft. 

The  velocity  of  the  plug  can  be  estimated  if 
values  of  M and  N are  known.  The  values  of  M and  N are 
not  well  known  for  delta  sediments.  General  estimates 
for  a Newtonian  viscosity  (i.e.,  N • 1)  have  been  set  at 
15  to  180  Ib-sec/sq  ft.  The  viscosity  of  subaer.'al 
debris  flows  has  been  estimated  for  a Bingham-Newtonian 
viscous  model  at  80,000  to  450,000  times  that  of  water. 


or  viscosities  3 to  15  Ib-sec/sq  ft. 


Estimates , 


however,  of  non-Newtonian  viscoelastic  properties  of 
clays  indicate  N of  about  0.07,  which  indicates  strong 
nonlinear  behavior.  ** 

As  an  illustration  of  the  method  for  determining 
the  flow  velocity  profile,  consider  the  case  for  N <*  1 
and  M • 10  Ib-sec/sq  ft.  The  Newtonian  model  for  the 
plug  velocity,  equation  (9),  becomes 


'8 


■ i (5^)(«o  - \)‘ 


and  the  veolcity  profile,  equation  (10), 


the  sides  of  the  gullies  could  be  expected.  The  thick- 
nesses of  the  plug  and  the  shear  zone  are  functions  of 
the  bottom  slope  and  the  pore  pressures  within  the 
sediments.  As  these  change  during  downslope  movement, 
the  thickness  of  the  layers  could  change.  Apparently  in 


the  toe  area  the  sheared  layer  decreases  in  thickness  as 
the  flow  decelerates,  so  that  when  stopped  the  toe 
thickness  reflects  the  plug  thickness. 


For  the  assumed  chute  characteristics  with  G • .99  , 
then  U “ 4.2  ft/sec:  and  for  G “ 1.0  , then  U = 6.S 

P . . P . 

ft/sec.  The  velocity  profile  is  shown  for  G “ .99  in 
Fig.  6.  The  flow  velocity  is  inversely  proportional  to 
the  viscosity  M,  so  that  doubling  M reduces  U to  one- 
half.  ^ 

Flow  of  the  sediments  in  the  toe  area.  The  toe 
areas  provide  the  first  opportunity  to  check  the  results 
of  the  theory  with  actual  data.  This  data  concerns  the 
shape  of  the  toe  noses  using  the  theory  of  Bruchl  and 
Scheidegger*"*  for  comparison.  The  toe  noses  display  a 
curving  edge  that  is  shown  in  Fig.  3.  The  first  example 
is  the  mudflow  toe  near  South  Pass  area.  A comparison 
between  the  observed  toe  shape  and  the  best  fit  A > 0.3 
to  the  theoretical  profile  is  shown  in  Fig.  7.  The 
values  used  in  the  example  were  h^  “ 60  ft,  y “ 90  Ib/cu 

ft,  P “ 0.021,  and  1))  • 20®.  The  best  fit  of  the  curve 

was  e • 0.022,  which  implies  C/C^  • 3.8  and  “ 119 

Ib/sq  ft.  This  value  of  e is  consistent  with  the  value 

at  failure,  and  the  C value  is  reasonable.  The  calcu- 

lated  ratio  of  C/C^  implies  C of  452  Ib/sq  ft,  a value 

that  is  higher  than  expected. 

A second  example  was  considered  from  the  Southwest 
Pass  area.  Selected  points  are  shown  in  Fig.  7.  The 
data  is  for  e ■ 0.02  and  h,  ■ 40  ft.  The  valuer 

D 

indicate  ” 72  Ib/sq  ft  and  a C/C^  ratio  of  2.8,  or  C ■ 
202  Ib/sq  ft.  The  best  curve  would  be  about  A * 0.2. 


Discussion 

The  results  of  this  study  seem  to  support  the 
application  of  a simple  rheological  model  to  the  mass 
movement  of  sediments  on  the  Mississippi  delta  front. 
The  rheological  model  used  combines  a static  Coulomb 
friction  stress  with  a viscous  stress,  much  in  the 
manner  used  to  describe  subaerial  mudflows.  The  data 
that  is  available  (mud  nose  shape)  seems  to  indicate 
general  agreement  with  theoretical  predictions.  Where 
data  is  not  available,  the  predictions  using  the  theory 
(velocity,  accelerations,  pore  pressures,  etc.)  seem  to 
be  reasonable,  or  at  least  correct  to  an  order-of- 
magnitude. 


As  encouraging  as  the  results  may  be,  they  must  be 
evaluated  in  light  of  the  present  general  lack  of  knowl- 
edge concerning  deformation  dynamics.  The  rheological 
model  used  is  not  the  only  one  that  could  have  been 
applied,  and  through  back  calculation  ("fudging"),  be 
forced  to  produce  reasonable  results.  Even  with  the 
model  used,  slight  changes  in  parameters  (i.e.,  viscos- 
ity M or  power  N)  that  are  poorly  known  would  greatly 
change  the  predictions,  making  them  at  times  indefen- 
sible. Another  limitation  to  these  results  stems  from 
the  assumptions  employed  in  applying  the  model,  i.e., 
steady  states,  constant  pore  pressure,  constant  cohe- 
sion, etc.  When  applying  this  approach  to  a specific 
site,  probably  any  specific  site,  these  assumptions  will 
almost  certainly  be  violated.  Finally,  the  data  used  in 
this  study  have  been  taken  from  several  sources  or  are 
outright  estimations  (values  for  M and  N) , so  that  a 
verification  of  the  predictions  cannot  be  claimed. 

These  results  indicate  that  further  studies  of 
submarine  mass  movements  should  consider  the  coordina- 
tion of  descriptive  studies  of  feature  geometries  (slump 
areas,  gullies,  toe  areas),  selection  and  placement  of 
in  situ  instruments  (spatially  and  with  depth),  and 
determination  of  geotechnical  (rheological)  properties 
of  the  sediments.  One  of  the  most  interesting  and 
potentially  important  prospects  for  future  study  is  that 
mass  movements,  as  exemplified  by  an  elongate  slide,  may 
be  a sequence  of  interrelated  events  separated  in  space 
by  thousands  of  feet,  extending  from  shallow  to  deep 
water  and  separated  in  time  by  hours,  if  not  days.  This 
possibility  means  that  the  value  of  making  measurements 
at  a single  site  to  explain  or  predict  sediment  movement 
at  that  same  site  may  be  quite  limited. 
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Table  1.  Pore  Pressure  Ratios  G at  Failure 
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F1g.  1 - Schematic  diagram  of  delta  front  Instability  features.  (From  Roberts  et  al 


ELONGATE  CHUTE 


MULTIPLE 

CROWN 

CRACKS 


Fig.  2 - Schematic  diagram  of  an  elongate  slide.  (From  Prior  and  Coleman  ) 


Fig.  3 •*  Sparker  record  showing  depositlonal  toe  and  the  rounded  frontal  scarp 
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Fig.  4 - Illustration  of  coordinate  axis  and 
geometry  of  soil  model . 


Fig.  5 - Theoretical  profiles  of  earth  flow  tongues  for 
various  values  of  the  parameter  A and  In  both  active  agd 
passive  Rankine  states.  (From  Bruchl  and  Scheldegger 


Fig.  6 - Example  of  vertical  velocity  profile  In  the  gully 
areas.  Calculations  based  upon  a viscosity  of  10  Ib-sec/sq  ft. 


Fig.  7 - Comparison  between  measured  mud  nose  shapes 
(dots  and  thick  line)  and  theoretical  profiles  from  Fig.  5. 


■Unclassified 

DOCUMENT  CONTROL  DATA  • R & D 

(S^<uritY  ftm»aHicmtion  oi  flff#,  body  ot  •batfct  mnd  inditing  mnnotmtion  ntuut  be  onffd  whmn  th»  oyf  >/  rmport  !•  clm$»lllmd} 

t.  ORICINATING  ACTIVITY  (CorporAlA  AuChoO  2A.  report  security  Ci.  ASSt  P 1 C A TION 

Coastal  Studies  Institute  ^ Unclassified 

Louisiana  State  University  26.  group 

Baton  Rouge,  Louisiana  70803  Unclassified 

3.  REPORT  TITLE 

EXPLANATION  OF  SUBMARINE  LANDSLIDE  MORPHOLOGY  BY  STABILITY  ANALYSIS  AND 
RHEOLOGICAL  MODELS  ^ 


4.  DESCRIPTIVE  NOTES  fType  o/ rApoff  •nct.IncluAivA 
5-  Au  THORISI  fFirsf  middtm  initiml,  tm0t  nmm^) 

Joseph  N.  Suhayda  and  David  B.  Prior 


• • REPORT  DATE 

May  1978 

•a.  contract  or  grant  no. 

N00014-75-C-0192 

b.  PROJEC  T NO. 

NR  388  002 


l7A.  TOTAL  NO.  OF  PAGES  176.  NO.  OF  REFS 


l»A.  ORIGlNATOR*S  REPORT  NUMBER(S) 


Technical  Report  No.  264 

•6.  OTHER  REPORT  NO(S)  (Any  othmr  numb^rm  thmt  may  ba  aaaignad 
thia  raport) 


10.  DISTRIBUTION  STATEMENT 

Approved  for  public  release;  distribution 

unlimited. 

II.  SUPPLEMENTARY  NOTES 

Reprint  from 

10th  Annual  OTC  in  Houston,  Tex.,  May 

1978,  pp.  1075-1082 

hXaosXract 

12.  SPONSORING  MILITARY  ACTIVITY 

Geography  Programs 

Office  of  Naval  Research 

Arlington,  Virginia  22217 

A theoretical  study  of  mass  movement  on  the  Mississippi  River  delta  front  has 
been  made  using  recently  acquired  field  data  and  a simple  rheological  model.  Recent 
measurements  of  sediment  properties  Include  cohesion,  bulk  density,  pore  pressure, 
and  internal  friction  angle.  Also,  the  geometry  of  a typical  type  of  instability 
feature,  an  elongate  slide,  is  examined  in  detail.  A rheological  model  describing  a 
Coulomb-viscous  plastic  in  effective  stress  terms  is  proposed  to  describe  certain 
features  of  the  mass  movement  process. 

The  model  is  used  to  derive  equations  defining  the  initial  failure  of  the  slopes, 
the  mass  movement  thickness  and  accelerations,  the  velocity  of  flow  in  the  gullies, 
and  the  shape  of  the  mud  nose  scarp  in  the  toe  area.  The  model  indicates  the 
Importance  of  pore  pressure  in  controlling  sediment  dynamics. 

Results  of  the  model  Compaq "fave^ably  with  the  observed  shapes  of  mud  noses. 
Initial  failures  on  low  slopes  (v 0.5°)^ are  attributed  to  pore  pressures  approaching 
geostatlc  values.  Flow  velociti&BTSffe  calculated  to  be  several  feet  per  second, 
based  upon  estimated  sediment  viscosities  and  pore  pressures  during  flow. 


/\6TT' 


DD .1473 

S/N  0101 -807-681  1 


(PAGE  1) 


Unclassified 
Security  Classification 


